Anisotropy of Solar Wind Turbulence between Ion and Electron Scales 
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The anisotropy of turbulence in the fast solar wind, between the ion and electron gyroscales, is directly 
observed using a multispacecraft analysis technique. Second order structure functions are calculated at different 
angles to the local magnetic field, for magnetic fluctuations both perpendicular and parallel to the mean field. 
In both components, the structure function value at large angles to the field S± is greater than at small angles 
Su: in the perpendicular component S±/S\\ — 5 ± 1 and in the parallel component S±/S\\ > 3, implying 
spatially anisotropic fluctuations, k± > ku . The spectral index of the perpendicular component is —2.6 at large 
angles and —3 at small angles, in broad agreement with critically balanced whistler and kinetic Alfven wave 
predictions. For the parallel component, however, it is shallower than —1.9, which is considerably less steep 
than predicted for a kinetic Alfven wave cascade. 
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Introduction. — Solar wind turbulence has been studied for 
many decades (e.g., and references therein) but a number 
of fundamental aspects of it remain poorly understood. This 
Letter will address one of these, the nature of the turbulent 
fluctuations at small scales, using a recently developed multi- 
spacecraft analysis technique. 

Turbulence is usually modeled as a cascade of energy, with 
injection at large scales and dissipation at small scales. In the 
solar wind, the injected energy is thought to originate from 
the observed large scale Alfven waves |2]. For scales between 
the effective outer scale and the ion gyroradius, termed the in- 
ertial range, a cascade of Alfvenic fluctuations H 3|— LTfl i s often 
invoked to explain the observed power spectra lla- ll ill . One 
aspect of recent investigation in the solar wind inertial range, 
relevant to this study, is anisotropy with respect to the mag- 
netic field. It has been shown that both power and searings 
vary with respect to the local magnetic field direction J |12l - ll6ll . 
in a way consistent with critical balance theories JUS]- 

At smaller scales, close to the ion gyroradius, the magnetic 
field power spectrum steepens (e.g., 11171 llal - ). This is com- 



anisotropy in the magnetic field components parallel and per- 
pendicular to the field are calculated. This provides a direct 
test of existing theories and a guide for new ones. 

Data set. — We use an interval of data from the Cluster mis- 



monly termed the dissipation range, although is sometimes 
called the dispersion range (e.g., B19I0 . and is where kinetic 
effects become important. Recent measurements of the mag- 
netic field spectral index in this range are between —2.3 and 
-2.8 JUHl |, although larger variation was seen in an ear- 
lier survey [18]. A further steepening in the spectrum near 
the electron gyroscale has also been observed J21, 23], In 
this study, we investigate between the ion and electron scales. 
Two popular suggestions for the types of fluctuations in this 
range are kinetic Alfven waves (KAWs) |0 on m " 
and whistler waves |T^, 25]. It has been suggested 17 
that, like some inertial range theories |5l|g]> the fluctuations 
are critically balanced, which would imply a spectral index 
of —7/3 in the perpendicular direction and —5 in the parallel 
direction. 

In this Letter, the first multispacecraft structure function 
measurements in the solar wind at scales below the ion gy- 
roscale are presented. The variance, power, and spectral index 
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sion B27I1 . in which the four spacecraft are in the fast solar 
wind with a separation ^100 km. The interval parameters are 
given in Table Hand are from the FGM JH, CIS & and 
PEACE 13011 instruments. No effects of Earth's foreshock are 
present, and the interval lies in the stable region of the param- 
eter space for pressure anisotropy instabilities (e.g., 113 ill ). 

For analyzing the fluctuations between ion and electron 
scales, high frequency data, > 1 Hz, is needed. In this study, 
a measurement of the local magnetic field direction is used, 
requiring data valid at both low and high frequencies. The 
STAFF instrument l32ll has a high frequency search coil mag- 
netometer, which provides a time series valid in the approxi- 
mate range 0.6-10 Hz. We combine this with the FGM data, 
which is valid up to «1 Hz in the solar wind. 

The combining procedure is based on the method in Ap- 
pendix A of Ref. [330- First, the high frequency (22 Hz) FGM 
data are interpolated onto the times of the STAFF data (25 
Hz). A wavelet transform is then applied to both time series 
to obtain two sets of wavelet coefficients for each field com- 
ponent. The wavelet coefficients corresponding to the STAFF 
data above 1 Hz are used to generate a high frequency time 
series, and those corresponding to the FGM data below 1 Hz 
are used to generate a low frequency time series. The two time 
series are then added, resulting in the combined signal. 

The power spectrum of the combined data is shown in Fig.Q] 
with the approximate ranges of FGM and STAFF marked. The 
noise floor of STAFF (from ground and in-flight tests 8331 ) is 
also shown. The break in the spectrum at «0.4 Hz is the ion 
scale spectral break point at the end of the inertial range, and 
is not due to the data merging. For this interval, the isotropic 
spectral index for the range of scales we study is — 2.88±0.01. 

Method. — A multispacecraft method is used in which data 
from the four Cluster spacecraft are combined to produce sec- 
ond order structure functions in different directions to the lo- 
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TABLE I: Interval parameters (Vsw, solar wind speed; n, number density; va, Alfven speed; T, temperature; p, gyroradius; d, inertial length; 
P, plasma beta) 
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cal magnetic field. It is based on the method of Ref. 11511 . One 
benefit of this technique is that a range of sampling angles can 
be covered simultaneously, enabling short intervals to be used, 
increasing the likelihood of statistical stationarity. 

Compared to the solar wind flow, the four spacecraft are 
approximately stationary and measure the magnetic field as 
the solar wind passes by. The measured variations, therefore, 
are due to both temporal and spatial variations in the plasma. 
Assuming that the temporal changes happen slowly compared 
to the flow, each time series can be converted into a spatial 
cut through the plasma (Taylor's hypothesis 0410 . Second 
order structure functions can then be calculated, defined as 
SB? (I) = (\Bi(r + 1) - B 4 (r)| 2 ), where B. L is the ith com- 
ponent of the magnetic field, 1 is the separation vector, and the 
angular brackets denote an ensemble average over positions r. 

It is important to consider the application of Taylor's hy- 
pothesis at small scales. In the inertial range, the solar wind 
speed is usually an order of magnitude larger than the Alfven 
speed and, therefore, Taylor's hypothesis is well satisfied. At 
smaller scales, the wave phase speed is larger than the Alfven 
speed l2lTl . It is still lower than the solar wind speed, 
however, so even if Taylor's hypothesis is less well satisfied 
than in the inertial range, it is not an unreasonable assump- 
tion. Measurements with an Alfvenic Taylor ratio (as defined 
in Ref. [15]) greater than 0.25 are discarded. 

Axisymmetry about the magnetic field is assumed so that 



• 10 
.10 2 

10° 
10" 
10" 

10- 

10" 

10"* 10" 




10"' 10"' 10" 10' 
Frequency (Hz) 



10 



FIG. 1: (color online). Magnetic field power spectrum of combined 
STAFF and FGM data from Cluster 3. Approximate instrument 
ranges are shown, as is the STAFF noise floor (green dashed line) 
and ion and electron gyroscales (blue dash-dot lines). Approximate 
frequencies corresponding to the range of scales used are marked 
(red dotted lines). 



1 can be split into parallel and perpendicular components, 
1 = (l\\,l±). There is mounting evidence that it is the lo- 
cal magnetic field that orders the fluctuations rather than a 
global field 11131 135143811 : i.e., the turbulence is anisotropic 
with respect to the field at the scale of each fluctuation rather 
than a much larger scale. Here, the local field is defined as 
B; oca ; = [B(r + 1) + B(r)]/2, and its direction is used to de- 
fine 2 1| and l±. The parallel and perpendicular components of 
B for each structure function pair are also defined with respect 

tO Blocal- 

The structure function values obtained from many pairs 
of magnetic field measurements from all four spacecraft are 
binned with respect to lu and l± and averaged. A minimum 
number of 200 values per bin is set to ensure reliable results 
and the binned data for each component are shown in Fig. |2] 
In both 8B\ and SB?, anisotropy can be seen: the structure 
function contours are elongated along the local field direction. 
Similar results have been seen at larger scales in inertial range 
solar wind correlation functions J39|], and in structure func- 
tions fromMHD iO and electron MHD JHHol] simulations. 

Instead of Cartesian coordinates, (l\\,lx), the data can be 
binned in polar coordinates, (1,9), where 9 is the angle to the 
local magnetic field (see Fig. [3]). Because of the low power, 
it is possible that the noise floor of STAFF has been reached 
for the small angle bins in SB%. This can be seen in Fig. [3] in 
which the lowest value structure function curves appear flatter 
than the others. Caution, therefore, is advised when interpret- 
ing these lowest power measurements. 

Variance anisotropy. — The ratio of power in the perpendic- 
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FIG. 2: (color online). Second order structure functions with respect 
to spatial separations parallel (lu) and perpendicular (l±) to the local 
magnetic field: 8B\ (left) and SB?, (right). 
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FIG. 3: (color online). Second order structure functions at different 
angles to the local magnetic field for the perpendicular component 
(upper set) and the parallel component (lower set). 



FIG. 4: (color online). Anisotropy of second order structure func- 
tions (upper) and spectral indices (lower). Various spectral index 
predictions are marked (dotted lines). 



ular component to the parallel component is sometimes re- 
ferred to as variance anisotropy, (e.g., j4llo . From Fig. [3] it 
can be seen that SB? is about 5% of SB\, which is smaller 
than average values of previous measurements [17, 4lh . This 
could be due to statistical variation, or due to the global, rather 
than local, mean field direction being used in those studies. 

The variance anisotropy for KAWs in electron reduced 
MHD 0| can be calculated, and for the parameters in Table 
Uthis prediction is 5B?/SB\ = 0.4, which is larger than the 
value observed here. Numerical solutions of linear kinetic the- 
ory, however, suggest smaller values of variance anisotropy 
for KAWs [42]. These values depend on f3, propagation an- 
gle, and wave number, but are in the range 0.01 to 0.2, which 
agrees with our result of 0.05. 

Power anisotropy. — Straight lines (in log-log space) are fit- 
ted to the data in Fig. [3] over the range 1 00^1-00 km. This is 
between ion and electron scales, i.e., between kpi = 1 and 
kp e = 1, where k — 2n/l. The interpolated values of the 
structure function at 200 km from these fits are given as a 
function of 9 in the upper panel of Fig. |4] The error bars are 
small, comparable to the data point size, and are the standard 
deviations of the points about the best fit lines. 

For both components, the structure function value 
("power") increases with 9. This is consistent with spatially 
anisotropic fluctuations, k± > kn, where k± and fey are char- 
acteristic parallel and perpendicular wave numbers fl43l1 . One 
measure of this anisotropy is the ratio of the largest angle bin 
value S± to the smallest S\\, which for SB\ is S±/S\\ = 5±1. 
This number is uncertain for SB?, due to the noise issues men- 
tioned above, but has a lower limit of 3. 



Some previous studies of the solar wind between ion and 
electron scales have measured anisotropy using the "slab plus 
2D" model Jl7ll4lll . The large slab fractions obtained are not 
generally in agreement with this study. In fact, our results are 
more consistent with other solar wind 11141 14411 and magne- 
tosheath j45ll measurements that demonstrate significant spa- 
tial anisotropy at scales smaller than the ion gyroscale. 

Spectral index anisotropy. — An important characteristic of 
turbulence is the spectral index —a of the power spectrum, 
E{k) ~ k~ a . Second order structure function scalings, i.e., 
gradients g, of the straight line fits to the data in Fig. [3] are 
related to the spectral index by a = g + 1 ll46ll . Using this re- 
lationship, the spectral index as a function of angle 9, for 5B\ 
and SB?, is shown in Fig. [4] The error bars are the standard 
errors on the best fit line gradients. 

For 8B\, the spectral index varies from around —2.6 at 
large angles to —3.2 at small angles. It should be noted that 
the steepest spectral index it is possible to measure with this 
method is —3 (e.g., |4(| 46]); for steeper spectra, the scal- 
ing seen by the two-point second order structure function is 
g = 2, since it is dominated by the smooth variation of the 
large scale field. At small angles, we observe a spectral in- 
dex of —3.2 ± 0.2, indicating that the spectrum in the parallel 
direction is k^ 3 or steeper. The predictions for a critically bal- 
anced whistler or KAW cascade are —7/3 in the perpendicu- 
lar direction and —5 in the parallel direction |7l|26[]. Although 
the spectral indices in Fig. Hare slightly steeper than the pre- 
diction at large values of 9, the steepening towards small 9 is 
suggestive of a critically balanced cascade. 



4 



The spectral index of SB?, varies from —1.9 at large an- 
gles to —1.5 at small angles. The small angle values may be 
affected by noise (as discussed previously), but the large an- 
gle ones appear not to be, and are significantly shallower than 
those of 5B\. This difference in gradient between the compo- 
nents can also be seen in Fig. [3] For a KAW cascade, SB?, is 
expected to scale in the same way as 5B\ Q2D - The difference 
observed here, therefore, may be indicating the presence of 
other modes or a different cascade mechanism. Another pos- 
sibility for the difference is instability generated fluctuations, 
although the measured parameters suggest the interval is not 
unstable to pressure anisotropy instabilities (e.g., 13111 '). 

Summary and conclusions. — The variance, power, and 
spectral index anisotropy are measured in the fast solar 
wind, between the ion and electron gyroscales. The vari- 
ance anisotropy is significant, with SB?^ being approximately 
5% of SBj^, Both magnetic field components display power 
anisotropy, implying spatially anisotropic fluctuations, k± > 
fcii. The spectral index of 8B\ steepens at small angles to the 
field, which is consistent with a critically balanced cascade of 
whistlers or KAWs. The spectral indices of <5£^are less con- 
sistent with the predictions, suggesting that the KAW picture 
|0] may be incomplete. 

Although we have looked for other data intervals, it is hard 
to find ones that satisfy the conditions required for this analy- 
sis, i.e., ~1 h long, away from Earth's foreshock, with small 
spacecraft separations and good angular coverage. A larger 
study is required to determine if the behavior noted here is 
typical for the solar wind. This may need to wait for a fu- 
ture mission due to the limitations of multispacecraft data cur- 
rently available. 
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